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Abstract 

Hyperfine induced 2s2p 3 Pq — > 2s 2 1 Sq transition rates in an external magnetic 
field for Be-like 47 Ti were calculated based on the multiconfiguration Dirac-Fock 
method. It was found that the transition probability is dependent on the magnetic 
quantum number Mp of the excited state, even in the weak field. The present 
investigation clarified that the difference of the hyperfine induced transition rate 
of Be-like Ti ions between experiment [Schippers et ah, Phys Rev Lett 98, (2007) 
033001(4)] and theory does not result from the influence of external magnetic field. 
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1 Introduction 



The hyperfine induced transition (HIT) rate of the 2s2p 3 Pq level for Be-like 47 Ti ions has 
been measured with high accuracy by means of resonant electron-ion recombination in the 
heavy-ion storage-ring TSR of the Max-Planck Institute for Nuclear Physics, Heidelberg, 
Germany [T]. However, the measured transition rate Ahit = 0.56(3) s" 1 differs from all 
present theoretical results Ahit ~ 0.67 s -1 [21 [3J, H] by about 20%. In the theoretical 
calculations the major part of the electron correlation, which always causes the dominant 
uncertainty, has been taken into account very elaborately. As a result, it is desirable to 
find out other reasons for the difference. 

In this letter, we focus on the influence of the magnetic field present in the heavy-ion 
storage-ring on the HIT rate. The HIT rate in an external magnetic field depends on the 
magnetic quantum number Mp of the excited state, even in a relatively weak field. This 
effect, combined with the non-statistical distribution of the magnetic sublevel population 
of the excited level, might lead to the difference in transition rate mentioned above. 

2 Theory 

In presence of the magnetic field, the Hamiltonian of an atom with non-zero nuclear spin 
I is 

H = Hf s + Hhf s + H m , (1) 

where Hf s is the relativistic fine-structure Hamiltonian that includes the Breit interac- 
tion. Hhf s is the hyperfine interaction Hamiltonian, which can be written as a multipole 
expansion 

%^T«.M( fc », (2) 

where and M( fc) are spherical tensor operator in electronic and nuclear space, respec- 
tively [5] . H m is the interaction Hamiltonian with the external homogeneous magnetic field 
B, 

H m = (N( 1} + AN (1) ) • B, (3) 

where are first-order tensor with the similar form of T/W, AN* 1 ' is the so called 
Schwinger QED correction [B]. 
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We choose the direction of the magnetic field as the ^-direction, and only Mp is a 
good quantum number. The wavefunction of the atomic system can thus be written as an 
expansion 

\TTIMp) = ^d TJF \TTIJFM F ). (4) 

TJF 

The total angular momentum F is coupled by the nuclear I and electronic J angular mo- 
mentum. The T and V are the other quantum numbers labeling the nuclear and electronic 
states, respectively. 

The coefficients drjF in Eq. fll]) are obtained through solving the eigenvalue equation 
using HFSZEEMAN package [7J 

Hd = Ed, (5) 
where H is the interaction matrix with elements 

Htjf,t>j>f> = (TTIJFM F \H fs + H hfs + H m \TT'IJ'F'M F ). (6) 

The readers are referred to Ref. 0, U\ for a detailed derivation of the different matrix 
elements . 

For the present problem, the wavefunction of the 3 Pq state can be written 
| u 2s2p 3 P / Mp") = d \2s2p 3 P I F(= I) M F ) + d s-,F'\ 2s2 P Sp i 1 F ' m f)- (7) 

5(=1,3);F' 

The quotation marks in the left-hand wave function emphasize the fact that the notation is 
just a label indicating the dominant character of the eigenvector. Remaining interactions 
between 2s2p 3 P and higher members of the Rydberg series can be neglected due to large 
energy separations and comparatively weak hyperfine couplings [S]. Furthermore, those 
perturbative states with different total angular momentum F can be neglected because of 
relatively weak magnetic interaction. As a result, Eq. ([7]) is simplified to 

| "2s2p 3 P / M F n ) = d \2s2p 3 P I F{= I) Mp) + ^ d s \2s2p S P X I F(= I) Mp). (8) 

5=1,3 

Similarly, the wavefunction of the ground state is approximatively written 

| "2s 2 / M F n ) = \2s 2 1 S I F(= I) M F ), (9) 
where all perturbative states were neglected for the same reasons as mentioned above. 
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The one-photon 2s2p 3 P — > 2s 2 1 Sq El transition becomes allowed via mixing with 
the perturbative states of 2s2p 3 P 1 and 2s2p 1 Pi (see Eq. OH])) induced by both the off- 
diagonal hyperfine interaction and the interaction with the magnetic field. The decay rate 
a(M F ) HIT from the excited state | "2s2p 3 P / M|") to the ground state | "2s 2 / Mf ") 
in s -1 is given by 

2.026 
"A 3 



a{M F ) HIT = 2 - U2bl6 x K"2s 2 / Mf"|P«|"2s2p 3 P I M F »)\ 2 , (10) 



Substitute Eq. OH]) and §S§ into above formula, then 

a(M F ) HIT = 2 ' 026i y 1018 £ I £ <W2F*(= J) + l v / 2F e (= + l 

9 5 



F»(= /) 1 P e (= J 9 {= 0) Fs(= I) I 

-M 9 F{=I) q M Fe[=I) ) { P e (= I) J e (= 1) 1 



(2s 2 1 S \\P w \\2s2p s P 1 )\ 2 . 

(11) 

Applying standard tensor algebra, the Eq. (II ip is further simplified to 



a(M' F ) H , T = 2026 ! j 3 A3 X WW (2I +1 )Y,\E d s ( _' M , \ ^J(2 S 21 S„||P( 1 »||2 S 2^P 1 )P 

(12) 

where A is the wavelength in A for the transition and (2s 2 1 S'o||P^||2s2p S P\) the reduced 
electronic transition matrix element in a.u.. 

From the Eq. (j!2p we can obtain the Einstein spontaneous emission transition proba- 
bility [9] 

A(M F ) HIT — HIT 
2.02613 x 10 18 



3A 3 



J2ds(2s 21 S \\P^\\2s2p s P 1 )\ 2 . (13) 

s 

It should be noticed that in present approximation of weak magnetic field, i.e., neglect- 
ing those perturbative states with different total angular quantum number P, the formula 
for the transition rate (see Eq. [TBI is similar to the one where the transition is induced 
by only hyperfine interaction [2] [3] . However, a significant difference exists in the mixing 
coefficients ds by virtue of incorporating the magnetic interaction into the Hamiltonian for 
the present work. 
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The electronic wavefunctions are computed using the GRASP2K program package [TO] . 
Here the wavefunction for a state labeled 7J is approximated by an expansion over jj- 
coupled configuration state functions (CSFs) 



In the multi-configuration self-consistent field (SCF) procedure both the radial parts of 
the orbitals and the expansion coefficients q are optimized to self-consistency. In the 
present work a Dirac-Coulomb Hamiltonian is used, and the nucleus is described by an 
extended Fermi charge distribution [TTJ. The multi-configuration SCF calculations are 
followed by relativistic CI calculations including Breit interaction and leading QED effects. 
In addition, a biorthogonal transformation technique introduced by Malmqvist [121 H3] is 
used to compute reduced transition matrix elements where the even and odd parity wave 
functions are built from independently optimized orbital sets. 

3 Results and discussion 

As a starting point SCF calculations were done for the configurations belonging to the 
even and odd complex of n = 2, respectively. Valence correlation was taken into account 
by including CSFs obtained by single (S) and double (D) excitations from the even and 
odd reference configurations to active sets of orbitals. The active sets were systemati- 
cally increased up to n ^ 5. The SCF calculations were followed by CI calculations in 
which core-valence and core-core correlations and the Breit interaction and QED effects 
were incorporated. Based on this correlation model, we calculated the hyperfine induced 
2s2p 3 P ->■ 2s 2 El transition rate for Be-like 47 Ti ions in absence of the magnetic field 
to Ah it — 0.66 s" 1 , where the experimental wavelength 346.99 A [H] was used to re-scaled 
the rateQ The value is in good agreement with the other theoretical results: Ahit — 0-67 
s" 1 by Cheng et al. [2] and A H it = 0.677 s _1 by Andersson et al. [3J. 

Recent theoretical calculations are all in disagreement with the experimental measure- 
ment A = 0.56(3) s _1 [1] by about 20%. It is hypothezised that the discrepancy results 
from the effect of magnetic field present in the storage ring. Actually, the magnetic field 

1 The nucleus of 47 Ti has the nuclear spin I = 5/2, nuclear dipole moment n = —0.78848 in and 
electric quadrupolc moment Q = 0.3 in barns [15j . 




(14) 
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effect has already been noticed and been discussed in previous experiment measuring the 
lifetime of the hyperfine state of metastable level 5c? 4 Di/ 2 for Xe + using the ion storage 
ring CRYRING at the Manne Siegbahn Laboratory (Stockholm) [16]. Returning to the 
present problem, experiment was conducted in the heavy-ion storage-ring TSR where the 
rigidity of the ion beam is given as B x p = 0.8533 T pQ, and the bending radius of the 
storage ring dipole magnets is p = 1.15m [17]. As a result, the magnetic field in the exper- 
iment has been 0.742 T. Considering the factual experimental environment, we calculated 
the hyperfine induced 2s2p 3 Pq — > 2s 2 1 Sq El transition rate of Be-like 47 Ti ion in the 
external magnetic field B=0.5 T, B=0.742 T and B=l T, respectively. With assistance 
of Eq. ( II 2p and Eq. (TT3"j) . we obtained the transition rate cl(M f )hit from the excited 
Zeeman state to the ground Zeeman state, the Einstein transition probability A(M f )hit 
of the excited state, and the corresponding lifetime r. Computational results are displayed 
in Table 1. As can be seen from this table, the transition rates A(M f )hit for each of 
the individual excited states "2s2p 3 P I M F n are obviously different because the mixing 
coefficients ds in Eq. ( TIB"]) depend on the magnetic quantum number Mp of the excited 
state. 

As can be found from Table 1, the lifetime of 3 Pq level is still not sensitive to the 
sublevel specific lifetimes, if the magnetic sublevels are populated statistically (the lifetimes 
t = J2 M e r(Mp)/(2I + 1) = 1.52s, 1.52s, 1.53s in the external magnetic field B=0.5T, 
0.742T and IT, respectively). In this case, the zero-field lifetime within the exponential 
error can be obtained, as made in Ref. [T], through only a fit of one exponential decay 
curve instead of 6 exponential decay curves with slightly different decay constants. To the 
contrary, in the experiment measuring the HIT rate of the 2s2p 3 P level of the Be-like 
Ti ion, the level concerned was produced through beam-foil excitation [18]. As we know, 
the cross sections with magnetic sublevels for ion-atom collision are different [T^J [20], and 
the magnetic sublevel population is in general not statistically distributed. Combining this 
fact with the M^-dependent HIT rate in an external field, the transition probability of 3 Pq 
level cannot be obtained by statistical average over all magnetic sublevel. However, we also 
noticed that an external magnetic field can lower the transition rate only for those magnetic 
sublevels with M F ^0. In other word, only if these specific magnetic sublevels with 
Mp ^ were populated, it is possible to explain or decrease the discrepancy between the 
measured and theoretical HIT rates for Be-like 47 Ti. In fact, such extreme orientation of the 
stored ions seems improbable by means of beam-foil excitation. Moreover, the experimental 
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heavy-ion storage-ring was only partly covered with dipole magnets (this fraction amounts 
to 13%) [17]. It further reduces the influence of magnetic field on the lifetime of level. 
Therefore, we still cannot clarify the disagreement between experimental measurement 
and theoretical calculations at present even though the influence of an external magnetic 
field was taken into account. 

4 Summary 

To sum up, we have calculated the hyperfine induced 2s2s 3 Po ~+ 2s 2 1 5'o El transition rate 
in an external magnetic field for each of the magnetic sub-hyperfine levels of 47 Ti 18+ ions 
based on the multiconfiguration Dirac-Fock method. It was found that the transition rate 
is dependent on the magnetic quantum number M F of the excited state, even in relatively 
weak magnetic fields. Considering the influence of an external magnetic field, we still did 
not explain the difference in the HIT rate of Be-like Ti ion between experiment and theory. 
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Table 1: Hyperfine induced 2s2p 3 Po — > 2s 2 1 S'o El transition rates in presence of magnetic 
field B=0.5 T, B=0.742 T and B=l T for Be-like 47 Ti ion. a represents the transition 
probability from the excited state "2s2p 3 P I to the ground state "2s 2 % / M|", 
A is the Einstein transition probability from the excited state "2s2p 3 P I My . r is the 
lifetime of excited state "2s2p 3 Pq I Mp" . The experimental wavelength (A) 346.99 A|14j 



was used in this calculations, where the influence of hyperfine interaction and magnetic 
field was neglected. 
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